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ABSTRACT 



Axisymmetric MHD simulations are used to investigate the Bondi-Hoyle accretion 
onto an isolated magnetized neutron star moving supersonically (with Mach number of 
3) through the interstellar medium. The star is assumed to have a dipole magnetic field 
aligned with its motion and a magnetosphcric radius i? m less then the accretion radius 
i?BH, so that the gravitational focusing is important. We find that the accretion rate 
to a magnetized star is smaller than that to a non-magnetized star for the parameters 
considered. Close to the star the accreting matter falls to the star's surface along the 
magnetic poles with a larger mass flow to the leeward pole of the star. In the case 
of a relatively large stellar magnetic field, the star's magnetic field is stretched in the 
direction of the matter flow outside of i? m (towards the windward side of the star) . For 
weaker magnetic fields we observed oscillations of the closed magnetosphere frontward 
and backward. These are accompanied by strong oscillations of the mass accretion rate 
which varies by factors ~ 3. Old slowly rotating neutron stars with no radio emission 
may be visible in the X-ray band due to accretion of interstellar matter. In general, the 
star's velocity, magnetic moment, and angular velocity vectors may all be in different 
directions so that the accretion luminosity will be modulated at the star's rotation 
rate. 



1 INTRODUCTION 



Accretion onto a supersonically moving star is an impor- 
tant astrophysical process which has been studied theoreti- 
cally and with simulations over many years (review by Edgar 
2004). This type of accretion has applications in many ar- 
eas. For example, the accretion by a compact star in a binary 
system where the gas supply is from a stellar wind or from 
a common envelope (Petterson 1978; Taam and Sandquist 
2000); and flows in regions of star formation (Bonnell et al. 
2001; Edgar and Clarke 2004). 

Supersonic accretion of a pressureless fluid onto a grav- 
itating center was investigated analytically by Hoyle and 
Lyttleton (1939) and an analytic solution was found by 
Bisnovatyi-Kogan et al. (1979). Later, Bondi and Hoyle 
(1944) developed a model including the fluid pressure. The 
Bondi-Hoyle model gives an explicit prediction for the ac- 
cretion rate. 

Supersonic accretion has been studied using numerical 
hydrodynamic simulations by Ruffert (Ruffert, 1994a; Ruf- 
fert and Arnett, 1994; Ruffert, 1994b). He studied the flow 
of gas past an accretor of varying sizes. For small accretors, 
the accretion rates obtained were in broad agreement with 
theoretical predictions. Accretion of matter to a nonmag- 
netized gravitating center was investigated by Bisnovatyi- 



Kogan & Pogorelov (1997). Pogorelov, Ohsugi, & Matsuda 
(2000) studied the Bondi-Hoyle accretion flow and its depen- 
dence on parameters but also for a nonmagnetized center. 

Here, we investigate Bondi-Hoyle accretion to magne- 
tized stars in application to the accretion of the interstellar 
medium (ISM) onto Isolated Old Neutron Stars (IONS)Q 
These slowly rotating neutron stars may be visible due ac- 
cretion of the interstellar medium. These stars may have dy- 
namically important magnetic fields. We use axisymmetric 
magnetohydrodynamics simulations to model the accretion 
onto the star. Earlier, we investigated cases of supersonic 
motion of strongly magnetized stars through the ISM in cases 
of slow rotation (Toropina at all 2001; hereafter T01) and 
fast rotation in the propeller regime (Toropina, Romanova, 
Lovelace 2006; hereafter T06). In these cases the magneto- 
spheric radius is greater then accretion radius i? m ^ Rbh 
and gravitational focusing is not important. Both cases cor- 
respond to magnetars with magnetic field B > 10 14 G. 

The present simulations are analogous to those in T01 



1 Isolated neutron stars have several stages in their evolution. 
Initially they rotate rapidly and may be observed as a radio pul- 
sars. As the star spins down and its period increases above several 
seconds it is thought that the pulsar mechnism stops working. 
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Figure 1. Sketch of the Bondi-Hoyle accretion onto a supcrsonically moving magnetized star. The thick black line is the shock wave 
and the red lines are the field lines of the star. 



and T06. However, here the star has a more typical neutron 
star magnetic field, B ~ 10 10 — 10 12 G. The magnetospheric 
radius is less then the accretion radius and gravitational 
focusing is important. A star gathers matter from the region 
inside the accretion radius. 

Sections 2 and 3 describe the physical situation and 
simulation model. In Sec. 4 we discuss the results of our 
simulations. In Sec. 5 we apply our results to more realistic 
neutron stars. Conclusions are given in Sec. 6. 



2 PHYSICAL MODEL 

A non-magnetized star moving through the ISM captures 
matter gravitationally from the accretion or Bondi-Hoyle 
radius 



Rbh = 



2GM 



3.8 x 10 —= — cm 



(1) 



where vioo = «/(100 km/s) is the normalized velocity of 
the star, c s the sound speed of the undisturbed ISM, and 
Mi. 4 = M / (IAMq) is the normalized mass of the star. The 
mass accretion rate at high Mach numbers M = v/c s >> 1 
was derived by Hoyle and Lyttleton (1939), 



Mhl =4tt(GM) 2 4 



7.3 x 10 8 -5— Mf. 4 S , (2) 



where p is the mass-density of the ISM and n — n/1 cm~ 3 
is the normalized number density. 

The Hoyle-Lyttleton model neglects the fluid pressure 
as mentioned. Later, Bondi and Hoyle (1944) extended the 
analysis to include the pressure. A more general formula for 
arbitrary Mach numbers M — v/c 3 was proposed by Bondi 
(1952), 

9 



¥ B h = irRlnpv = 47ra(GM) 2 



(v 2 + c'§) s / 2 ' ^ 

where the coefficient a is of the order of unity. Bondi sug- 
gested a = 1/2, but here we assume a = 1 in agreement with 
equation (2). Values of the a parameter for different specific 
heat ratios are determined by the simulations of Pogorelov, 
Ohsugi, & Matsuda (2000). A sketch of the geometry for the 
Bondi-Hoyle analysis is shown on Figure 1. 

If the fast moving star has a significant magnetic field, 
then the situation is more complex. In the limit of a strong 



stellar magnetic field, the ram pressure of ISM can be bal- 
anced by the magnetic pressure of the star's magnetosphere 
at a radius i? m which is larger than Rbh- In this limit 
the ISM interacts directly with the star's "large" magne- 
tosphere. This case investigated earlier by T01 is relevant to 
stars with very strong magnetic fields (e.g., magnetars) or 
to stars with very high velocity. 

In the opposite limit where i? m < 7?bh the situation 
is different. First, matter is captured by the gravitational 
field of the star as in Bondi-Hoyle accretion onto a non- 
magnetized star and it accumulates. If the magnetosphere 
of the star is very small, R m <C -Rbh, and the gravitation- 
ally attracted matter accretes to the star in a spherically 
symmetric flow, then 



( BtRl \ 
\MV2GM J 



2/7 



(4) 



(e.g., Lipunov 1992), where B t is the magnetic field at the 
surface of the star of radius R, and M is the accretion rate. 

If the magnetized star accretes matter with the same 
rate as a non- magnetized star, M = A/bh, then the magne- 
tospheric radius is 



R n 



4.7 x 10 



R 4/7 R 12/7 6/7 
9 £>!() fMj ^100 



(5) 



from equation (4) with B w = _B„/10 10 G and R 6 = 
-R*/10 6 cm. For the adopted reference parameters, _R m is 
roughly equal to -Rbh/10 3 . Our simulations discussed be- 
low show that the accretion rate to the star is in general 
less than A/bh. 

Our earlier studies of accretion onto magnetized stars 
(Toropina at al. 2001; Romanova at al. 2003; Toropina at 
al. 2003) showed that the magnetosphere acts as an obstacle 
for the flow and that the star accretes at a lower rate than 
a non- magnetized star for the same v, c s , and M. In all 
those simulations we considered strongly magnetized stars 
including magnetars. 

In this work we consider stars with magnetic fields 
B ~ 10 10 — 10 12 G, and investigate its interaction with the 
interstellar medium. We neglect rotation of the star. 
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3 NUMERICAL MODEL 

We investigate the interaction of fast moving magnetized 
star with the ISM using an axisymmetric, resistive MHD 
code. The code incorporates the methods of local itera- 
tions and flux-corrected-transport developed by Zhukov, 
Zabrodin, & Feodoritova (1993). This code was used in our 
earlier simulations of Bondi accretion onto a magnetized star 
(Toropin at. al. 1999, hereafter T99), propagation of a mag- 
netized star through the ISM (T01), spherical accretion to a 
star in the "propeller" regime (Romanova at al. 2003; here- 
after R03) and spinning-down of moving magnetars in the 
propeller regime (T06). 

The flow is described by the resistive MHD equations: 



dp 
dt 



+ V-(p v) 







dv 1 
p— + p(v • V)v = -Vp + - J x B + F 9 



^=Vx(vxB) + -l-V 2 B 



dt 



4-kct 



+ V ■ {pew) = -pV-v + — 
at (j 



(6) 



We assume axisymmetry (d/d(j> = 0), but calculate all 
three components of velocity and magnetic field v and B. 
The equation of state is taken to be that for an ideal gas, 
p — (7 — l)pe, where 7 = 5/3 is the specific heat ratio and 
e is the specific internal energy of the gas. The equations 
incorporate Ohm's law J = <r(E + v x B/c), where a is 
the electric conductivity. The associated magnetic diffusiv- 
ity, rjm = c 2 /(47rer), is considered to be a constant within 
the computational region. In equation (2) the gravitational 
force, F 9 = — GMpR/R 3 , is due to the star because the self- 
gravity of the accreting gas is negligible. 

We use a cylindrical, inertial coordinate system (r, cf>, z) 
with the z— axis parallel to the star's dipole moment fi and 
rotation axis CI. The vector potential A is calculated so that 
automatically V • B = at all times. We neglect rotation 
of the star, ft — 0. The intrinsic magnetic field of the star 
is taken to be an aligned dipole, with vector potential A = 
H x R/R :i . 

We measure length in units of the Bondi radius Rb = 
GM/Cg — A4 2 Rbh/2, with c B the sound speed at infinity, 
the density in units of the density of the interstellar medium 
at infinity, p, and the magnetic field strength in units of Bo 
which is the field at the pole of the numerical star. Here 
we use the term numerical star to distinguish the model 
star (inner boundary) from the actual neutron star. The 
magnetic moment is measured in units of po = Bo-Rb/2. 

After reduction to dimensionless form, the MHD equa- 
tions involve the dimensionless parameters: 



= 



87rp() 



RbVq Re n 



(7) 



where ?) m is the dimensionless magnetic diffusivity, _Re m is 
the magnetic Reynolds number, and j3 is the ratio of the 
ISM pressure at infinity po and the magnetic pressure at the 
poles of the numerical star. 

Simulations were done in a cylindrical region (Zmin < 



z ^ Z max ,0 ^ r ^ Rmax). The numerical star was repre- 
sented by a small cylindrical box with dimensions R* « 
-Rmax and |Z„ I << Z max . A uniform (r, z) grid with N R x N z 
cells was used. 

Initially, the density p(r, z) and the velocity v(r, z) are 
taken to be constant in the region p(r, z) — p and v = v* , 
v$ — 0. Also, initially the vector potential A was taken to 
be that of a dipole so that B$ = 0. The vector potential was 
fixed inside the numerical star and at its surface during the 
simulations. 

The outer boundaries of the computational region were 
treated as follows. Supersonic inflow with Mach number M 
was specified at the upstream boundary (z — Z m in, ^ r ^ 
iJmax). The variables (p, v r , v z , e) are fixed with p and s 
corresponding to the ISM values, v z = 0, v r — Vao = A4c 3 . 
The inflowing matter is unmagnetized with B = 0. At the 
downstream boundary (z = Z max , < r < Rmax) and at the 
cylindrical boundary (Z m i n ^ z ^ X m „,r = i? max , a free 
boundary condition was applied, d/dn = 0. The boundary 
conditions are described in greater detail in T01 and T06. 

The size of the computational region for most of the 
simulations was J? m ax = 0.27, Z m in = —-Rmax = —0.27, with 
Z mal = 2Rmax = 0.54 in units of the Bondi radius. The 
theoretically derived accretion radius (Bondi-Hoyle radius) 
Rbh = 0.2. Therefore the size of the computational region 
is sufficient for our simulations. 

The grid N R x N z was 1297 x 433 in most of cases. The 
radius of the numerical star was R» = 0.0025 in all cases. 



4 RESULTS 

We consider Bondi-Hoyle accretion onto a star with a small 
magnetosphere. Simulations were done for a number of val- 
ues of the star's magnetic moments p. The numerical star 
(inner boundary) is typically much larger than the true ra- 
dius of the neutron star. Thus the actual star is "hidden" 
inside the numerical star. The Mach number for all of the 
present simulations is M. = 3. Also, f) = 10~ 6 , which corre- 
sponds to the magnetic moment p = 10 30 Gcm 3 . The mag- 
netic diffusivity is r/ m = 10 -5 . We measure time in units of 
to = {Zmax — Zmin )/v 00 = 0.81Rb/«oo, which is the cross- 
ing time of the computational region with the star's velocity, 



4.1 Classical Bondi-Hoyle Accretion 

To test our model we performed hydrodynamic simulations 
of the Bondi-Hoyle accretion to a non-magnetized star for 
Mach number M = 3. Dimensionless accretion radius for 
this parameter Rbh = 0.2. The star collects matter from 
inside this radius. We verified that the flow is close to that 
described earlier in T01 where matter forms a conical shock 
wave. 

Figure 2 shows the main features of the flow at a late 
time t — 7.0to when the flow is stationary. We calculated 
the accretion rate M to the numerical star and got a value 
M w 0.7Mbh for numerical star R, = 0.0025 were M B h is 
the Bondi-Hoyle accretion rate of equation (3) with a = 1. 
In our earlier calculation we got a slightly smaller value, 
M w 0.5M B h for a numerical star R» = 0.02 (T01). This 
difference can be explained by the higher resolution of the 
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Figure 3. Matter flow around a weakly magnetized star moving through the interstellar medium with Mach number Ai = 3 at a late 
time t = 4.4to. The background represents the logarithm of density and the solid lines are streamlines. The length of the arrows is 
proportional to the poloidal velocity 



present calculations. And this difference agrees with Ruf- 
fert's results on the dependence of M on numerical star size 
for the sizes used (Ruffert 1994). 

4.2 Bondi-Hoyle accretion to a star with a small 
magnetosphere 

We performed simulations for a number of values of the mag- 
netic moment. Initially we chose fi = 10 30 Gcm 3 (dimension- 
less moment /i = 1.6 x 10 -8 ), which corresponds to a star's 
surface magnetic field B = 10 12 G. We investigated this case 
in detail. 

Figure 3 shows the main features of the flow. Because of 
the magnetospheric radius R m << Rbh, the gravitational 



focusing is important and the star gathers matter from the 
region inside the accretion radius -Rbh ~ 0.18, which cor- 
responds approximately to the theoretically derived accre- 
tion radius (Bondi-Hoyle radius) -Rbh = 0.2. The magne- 
tospheric radius R m — 0.02 is about 10 times smaller than 
the gravitational radius Rbh- At the same time the average 
magnetospheric radius R m is about 3 times larger than the 
numerical star radius R» = 0.0025. The star's magnetic field 
at distances r > R m is significantly affected by the matter 
flow. An excess of matter flows from the backside of the star 
(return flow) and consequently the magnetic field lines are 
pushed towards the "front" of the star forming "ear" -like 
shapes. The "ears" oscillate in the z— direction. The conical 
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Figure 4. A central region around a weakly magnetized star moving through the interstellar medium with Mach number Ai = 3 at a 
late time t = AAto. Zoomed- in view of the stellar magnetic field which is deformed by the flow. The thick black line surrounding the star 
is the magnetospheric radius r m . The gray circle is the numerical star. 



shock wave is close to that found in the hydrodynamic case 
which is expected since the Mach numbers are the same. 

Figure 4 shows zoomed-in view of the stellar magnetic 
which is deformed by the flow at a late time t = 4.4£o- The 
thick black line surrounding the star is the magnetospheric 
radius r m where pv 2 = B 2 /87r. 

Figure 5 shows the dependences of the density, veloc- 
ity and temperature along the z— axis. There is a normal 
shock wave at z « —0.04. The density increases as the star 
is approached; and it is zero in the region of the star where 
matter is accreted. Behind the star (z > 0) there is a strong 
density enhancement owing to the gravitational focusing. 
At t = 4.4to the density of matter behind the star is signifi- 
cantly larger than the density on the front side. 

The velocity v z decreases sharply in the shock wave to 
subsonic values, but later increases again in the polar col- 
umn. Behind the star the velocity is negative in the region, 
where accretion occurs and also increases again in the po- 
lar column to supersonic values. The plasma temperature T 
and pressure p increase strongly approaching the star along 
the ±z directions. This is due to the adiabatic compression 
of the matter. 

Figure 6 shows the variation of the plasma pressure p, 
magnetic pressure p mag and p = p/p mag along the z— axis. 
The dashed lines correspond to the initial distributions. One 
can see that the pressure p gradually increases along z— axis, 
has a maximum in polar columns, and then gradually de- 
creases. The magnetic pressure p m a g initially has a dipole 
shape, but it is significantly modified by the inflow of mat- 
ter. Outside of the magnetosphere the matter pressure dom- 
inates. 

4.3 Oscillations 

We performed a number of simulation runs with different di- 
mensionless magnetic moments ranging from fj, = 4.7 x 10 -9 
to /i = 1.2 x 10 -7 . For relatively small magnetic moments, 
/i < 2 x 10 -8 , we observed oscillations of the magnetic field 



configuration as shown in Figure 7 for our main case where 
/i = 1.6 x 10~ 8 . These oscillations are connected with the 
balance between back flow of the accreting matter and the 
tension in the distorted stellar magnetic field. Figure 7 il- 
lustrates time evolution of magnetic field lines. The field 
oscillations are accompanied by strong oscillations in the 
accretion rate to the star as shown in the top panel of Fig- 
ure 8 for our main case where fi = 1.6 x 10~ 8 . For larger 
values of \x the oscillations are much weaker as shown in the 
bottom panel of Figure 8. 

A Fourier analysis of the accretion rate oscillations for 
our main case gives the quasi-period of these oscillations 
AT « 0.2t = 0.16i?B/i>oo = 0.72i? BH /«oo. 

Note that in hydrodynamic simulations Ruffert (1994b, 
1996) observed a "flip-flop" instability of the accretion flow 
in 2D simulations and mass accretion rate fluctuation for 
smaller accretors in 3D simulations. In our simulations we 
see instability of the flow in case of weak or moderate mag- 
netic field. This instability disappears as the magnetic field 
increases. Further, this instability is absent in our purely 
hydrodynamic simulation with the same parameters. 



4.4 Accretion rate 

For our main case the mass accretion rate is approximately 
M fa 0. 45Mb h- This value is smaller than mass accretion 
rate in the hydrodynamic case with the same parameters 
M w 0.7 Mb h- Hence even a small stellar magnetic field 
acts to reduce the mass accretion rate. This is in agreement 
with the conclusions of our earlier work (T99, T03). 

Figure 9 shows the observed dependence of the accre- 
tion rate on the star's magnetic moment. It is approximately 
M ~ fi~ 0A . Note however that as fj, increases, the magne- 
tospheric radius R m can approach the accretion radius -Rbh 
and this can suppress the accretion. 
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Figure 5. Panel (a) shows the density, panel (b) the axial velocity and panel (c) the temperature variation along the z— axis for our 
main case at t = 4.4to- The dashed lines give the initial distributions. 




Figure 6. Panel (a) shows the matter pressure p, panel (b) the magnetic pressure p m ag, and panel (c) = p/pmag variation along the 
2— axis for our main case at t = 4.4Tb. The dashed lines give the initial distributions. 
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Figure 7. Matter flow around a weakly magnetized star moving with Mach number JVi = 3 at the times indicated on the plots. The 
background color represents the logarithm of density and the length of the arrows is proportional to the poloidal velocity. The solid lines 
are magnetic field lines. Distances are measured in units of the Bondi radius. 



5 ASTROPHYSICAL PARAMETERS 

Let us consider an isolated old neutron star (IONS) mov- 
ing through the interstellar medium. The star is magnetized 
and has a mass M = 1.4 M©. The density of the inter- 
stellar matter is assumed to be p^, = 1.7 x 10~ 24 g/cm 3 
(«„ = 1/ cm 3 ). The sound speed in the ISM is assumed to 
be Cs = 12km/s so that at a Mach number M = 3 the star 
moves with a velocity v = 36 km/s which corresponds to 
a relatively slow moving neutron star (Arzoumanian, Cher- 
noff, & Cordes 2002). 

The Bondi-Hoyle radius is Rbh = GN'L/v 2 « 1.3 x 
10 13 cm. The Bondi-Hoyle accretion rate is the Mbh ~ 5 x 
lO~ 17 M0/yr for a non-magnetized star. This gives an accre- 
tion energy release rate of L acc = GMM / R* « 6 x 10 29 erg/s 
for R„ = 10 6 cm and a non-rotating star. Both M and 
L acc are reduced for magnetized star as as shown in Fig- 
ure 10. With ft = StvPoc/Bq and — piSMct/j, we ob- 
tain the magnetic field at the surface of the numerical star 
Bo = {8-K/-yf3) 1/2 p 1/2 c s w 0.006G. 

Recall that the size of the numerical star is R, = 
0.0025-Rb ~ 3 x 10 11 cm in all simulations. Thus the mag- 
netic moment of the numerical star is p = B rl/2 w 10 30 
Gem 3 . Note that this is of the order of the magnetic moment 
of typical spin powered radio pulsars which have a surface 
field L?ns ~ 10 12 G and radius of the order of 10 6 cm. Thus 



we can suggest that an actual neutron star exists inside of 
our numerical star. 

We can estimate the radius of the magnetosphere using 
equation (5) as _R m w 6.2 x 10 1 cm. Thus R m << -Rbh- 



6 CONCLUSIONS 

We carried out axisymmetric MHD simulations of Bondi- 
Hoyle accretion onto a fast moving neutron star (M = 3) 
with an aligned dipole magnetic field for cases where the 
Bondi-Hoyle accretion radius -Rbh is much larger than the 
star's magnetospheric radius R m . From this we conclude 
that: 

1. Interstellar matter is captured by the star if it has an 
impact parameter less than the accretion radius -Rbh - It then 
falls onto the star mainly onto the magnetic poles of the star. 
More matter falls downstream of the star. The accretion rate 
to the star is smaller than that to a non-magnetized star. 

2. In case of a strong stellar magnetic field (Rm > Rbh) 
the field is stretched in the direction of the dominant flow 
which towards the front side of the star. 

In the case of a weak magnetic field (Rm < Rbh), 
the star's field oscillates in the front/back direction. These 
oscillations are accompanied by strong oscillations of the 
mass accretion rate. The time scales of the oscillations is 
~ 0.7-Rbh/u CO • 
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Figure 8. The top panel shows the Dependence of mass accretion rate M on time for our main case with fi = 1.6 X 10 — 8 and the bottom 
panel a case with a stronger magnetic field, fi = 1.2 X 10 — 7 The dashed lines give the mass accretion rate normalized in Bondi-Hoyle 
rate, while the solid lines give the integrated mass flux. Time is measured in the crossing time units, Az/i^oo. 




Figure 9. Dependence of the mass accretion rate M on magnetic moment fi. The straight line is a least-squares straight line fit to the 
points. 



3. The accretion rate decreases with increasing mag- 
netic moment approximately as M ~ /i -0 ' 4 . 

4. Old slowly rotating neutron stars with no radio emis- 
sion may be visible in the X-ray band due to accretion of 
interstellar matter. For a neutron star moving through the 
interstellar medium with Mach number of 3 we estimate 
the accretion luminosity as < 10 30 erg s _1 . In general, the 
star's velocity, magnetic moment, and angular velocity vec- 
tors may all be in different directions so that the accretion 
luminosity will be modulated at the star's rotation rate. 
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